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We investigate the critical dynamics of the n-component relaxational models C and D which 
incorporate the coupling of a nonconserved and conserved order parameter S, respectively, to the 
conserved energy density p, under nonequilibrium conditions by means of the dynamical renormal- 
ization group. Detailed balance violations can be implemented isotropically by allowing for different 
effective temperatures for the heat baths coupling to the slow modes. In the case of model D with 
conserved order parameter, the energy density fluctuations can be integrated out, leaving no trace of 
the nonequilibrium perturbations in the asymptotic regime. For model C with scalar order parame- 
ter, in equilibrium governed by strong dynamic scaling (zs = z p ), we find no genuine nonequilibrium 
fixed point either. The nonequilibrium critical dynamics of model C with n = 1 thus follows the 
behavior of other systems with nonconserved order parameter wherein detailed balance becomes ef- 
fectively restored at the phase transition. For n > 4, the energy density generally decouples from the 
order parameter. However, for n = 2 and n = 3, in the weak dynamic scaling regime (zs < z p ) entire 
lines of genuine nonequilibrium model C fixed points emerge to one-loop order, which are character- 
ized by continuously varying static and dynamic critical exponents. Similarly, the nonequilibrium 
model C with spatially anisotropic noise and n < 4 allows for continuously varying exponents, yet 
with strong dynamic scaling. Subjecting model D to anisotropic nonequilibrium perturbations leads 
to genuinely different critical behavior with softening only in subsectors of momentum space and 
correspondingly anisotropic scaling exponents. Similar to the two-temperature model B (randomly 
driven diffusive systems) the effective theory at criticality can be cast into an equilibrium model D 
dynamics, albeit incorporating long-range interactions of the uniaxial dipolar or ferroelastic type. 

PACS numbers: 64.60.Ak, 64.60.Ht, 05.40.-a, 05.70.Jk 
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I. INTRODUCTION 



Analytical studies of dynamic critical phenomena in 
the vicinity of a second-order phase transition usually 
rely on a coupled set of Langevin-type stochastic equa- 
tions of motion for the relevant slow variables, namely 
the order parameter and hydrodynamic modes associ- 
ated with conservation laws pj. Taking advantage of 
the separation of time scales induced by critical slowing 
down, all remaining microscopic degrees of freedom are 
reduced to additive Gaussian white noise terms in this 
description. In order to guarantee that the probability 
distribution for any configuration converges to the canon- 
ical Gibbs function V cq (T) = Z(T)- X exp{-H/k B T) at 
long times t — > oo (with the effective Hamiltonian Ti. 
usually taken to be the standard 4 model), the second 
moments of the stochastic forces must be related to the 
relaxation rates via Einstein relations. In addition, inte- 
grability conditions constrain the reversible force terms 
in the nonlinear Langevin equations quite severely, for 
the associated probability currents in the space of the 
slow variables must be divergence- free Q- These two 
requirements also ensure the validity of the equilibrium 
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fluctuation-dissipation theorem which relates the imagi- 
nary part of the dynamic susceptibilities with the corre- 
lation functions. As a consequence, the system's static 
behavior can be separated from its dynamic properties. 

In isotropic systems, there are normally two indepen- 
dent static critical exponents, e.g., v and r\ which re- 
spectively characterize the divergence of the correlation 
length upon approaching the transition, £ ~ |r| ~ v , where 
r oc T — T c , and govern the power-law decay of the two- 
point correlation function at T c , C(|x|) ~ |x|~( d ~ 2+I ') 
in d spatial dimensions. These become supplemented by 
dynamic exponents z that describe the critical slowing 
down for the relevant modes, with characteristic relax- 
ation times diverging as t c ~ \t\~ zv . At thermal equi- 
librium, the dynamic universality classes are well under- 
stood, and known to be distinguished by overall features 
of the dynamical system at hand. In addition to the order 
parameter symmetry which essentially dictates the static 
critical exponents, the determining factors are if the or- 
der parameter itself represents a conserved quantity or 
not, the absence or presence of additional conservation 
laws, and the form of the reversible mode couplings be- 
tween the generalized hydrodynamic variables, as again 
dictated by the symmetries of the problem 0. 

On the other hand, critical dynamics in systems far 
from thermal equilibrium is not subject to the stringent 
limitations imposed by the detailed balance constraints, 
and in fact cannot even always be adequately captured 
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through coarse-grained stochastic equations of motion 
0- Nevertheless, several important situations have been 
successfully modeled by means of the Langevin formal- 
ism, two prominent examples being driven Ising lattice 
gases, or more general driven diffusive systems and 
nonequilibrium interface growth models, such as cap- 
tured by the Kardar-Parisi-Zhang (KPZ) equation and 
its variants |5| . Yet in nonequilibrium circumstances one 
has to invoke heuristic and/or phenomenological argu- 
ments for choosing the mathematical form of the noise 
correlations. This must be done with appropriate care, 
however, since the structure of the stochastic forces may 
crucially impact the universal scaling behavior [(|. It 
is thus of vital importance to elucidate the influence of 
different forms of the assumed stochastic noise correla- 
tions on the long-distance and long-time properties of any 
nonequilibrium Langevin system under investigation. 

Naturally then, the following question arises for the 
Langevin models describing the equilibrium dynamical 
universality classes: What happens to their universal 
scaling behavior if the detailed balance conditions are 
violated ? Notice that since 'static' properties cannot be 
decoupled from the dynamics in nonequilibrium steady 
states, this may include novel static critical behavior 
as well in addition to perhaps modified values for the 
dynamic exponents z. The simplest dynamical model 
just entails a purely relaxational kinetics for a noncon- 
served order parameter with no coupling to other con- 
served quantities. This defines model A in the (alpha- 
betical) classification of Ref. p|. Yet the model A uni- 
versality class, as for example realized in the kinetic Ising 
model with Glauber spin flip dynamics, is known to be 
extremely robust against nonequilibrium perturbations 
0, ■ For the kinetic Ising model, this remains true even 
when the order parameter up-down symmetry is broken 
@. Consider the most straightforward situation where 
the order parameter symmetry remains preserved, but 
the Einstein relation is not satisfied. Since there is only 
a single stochastic equation of motion in this case, one 
can recover detailed balance through simple parameter 
rescaling which does not affect universal properties at 
the phase transition [T(ij ]. 

Remarkably, the situation is markedly different for the 
purely diffusive relaxational critical dynamics of model 
B with conserved order parameter (e.g., the kinetic Ising 
model with Kawasaki spin exchange processes) , but only 
when subject to spatially anisotropic noise, say with 
stronger noise correlations in the thus defined longitudi- 
nal as compared to the complementary transverse sector 
in momentum space. In this effective two-temperature 
or randomly driven model B, excitations in the trans- 
verse sector soften first, while the longitudinal directions 
remain noncritical IT2I Il3| . This induces inherent 
anisotropic scaling at the critical point, of the same form 
as those in driven lattice gases Q. Interestingly, the 
emerging long- wavelength dynamics in the critical regime 
can be recast into an equilibrium model B, albeit with an 
effective Hamiltonian that incorporates long-range inter- 



actions of the uniaxial dipolar or ferroelastic type. These 
reduce both the lower critical dimension, allowing long- 
range order already in one dimension, as well as the up- 
per critical dimension to d c = 4 — d\\, where d u denotes 
the dimension of the stiff longitudinal sector [111 [l2| . In 
Ref. ^| , the associated four independent critical expo- 
nents (is, r/, z, and the anisotropy exponent A) were com- 
puted to two-loop order in the e expansion (e = d c — d) 
by means of the dynamic renormalization group (RG), 
utilizing a path integral or dynamic field theory repre- 
sentation of the Langevin equation Q . 

The equilibrium dynamical models C and D (in the ter- 
minology of Ref. 0) still describe purely relaxational dy- 
namics for either a nonconserved (model C) or conserved 
(model D) n-component order parameter field S, which is 
however statically coupled to the conserved scalar energy 
density p [l^. Preserving the 0(n) order parameter sym- 
metry, the lowest-order coupling is cx p S 2 . As the energy 
density itself represents a noncritical variable entering the 
Hamiltonian only quadratically, it can be integrated out 
exactly in the partition function Z(T) and the generat- 
ing function for static correlations. This merely shifts the 
value of the fourth-order coupling u for the order param- 
eter fluctuations, whence one recovers the static critical 
exponents of the 0(n) model. 

The coupling to the scalar diffusive mode p may how- 
ever alter the dynamic critical behavior. To one-loop 
order, three distinct scaling regimes emerge for model 
C with nonconserved order parameter, depending on the 
component number n pL5L Il6l fl7| | : (a) for Ising symmetry 
(n = I), one finds 'strong' scaling, i.e., zs — z p = 2+a/v, 
where a denotes the specific heat critical exponent; (b) 
the interval 2 < n < 4, for which a > 0, is characterized 
by 'weak' scaling with zs = 2(1 + a/nu) < z p = 2 + a/v; 
(c) for n > 4, where a < 0, the Langevin equations for 
the S and p effectively decouple, leaving purely diffu- 
sive behavior for the conserved mode, z p = 2, and the 
model A dynamic critical exponent for the order param- 
eter, z s = 2 + cr/, with c = 61n| - l + 0(e = 4-d). To 
higher orders in perturbation theory, these three regimes 
essentially persist (yet there appear additional distinc- 
tions with respect to the corrections to the leading scal- 
ing laws), but their boundaries become functions of the 
spatial dimension d as well as of n [Til Il7| . For model D 
with conserved order parameter, the energy density al- 
ways fluctuates faster in the critical region, rendering a 
strong-scaling regime impossible. The order parameter 
dynamics is thus not affected by the additional conserva- 
tion law, and given by the model B dynamic critical expo- 
nent zs = 4 — ?7. For a > 0, one finds again z p = 2 + a/ v, 
whereas z p = 2 in the decoupled case when a < |15|. 

In this paper, we explore the effect of perturbations 
in the stochastic force correlators that violate the equi- 
librium conditions on the critical dynamics of the relax- 
ational models C and D. Specifically, we shall retain the 
0(n) order parameter symmetry, but introduce differ- 
ent noise correlation strengths for the critical fluctuations 
and the conserved energy density, respectively, amount- 



3 



ing to unequal effective heat bath temperatures T5 and 
T p . We shall employ the dynamic RG to one-loop order, 
and search for novel nonequilibrium fixed points of the 
ensuing RG flow equations. In addition, we will investi- 
gate spatially anisotropic detailed balance violations. 

The critical dynamics at structural phase transitions 
and of anisotropic antiferromagnets are usually listed as 
possible realizations of the model C universality class . 
In the latter case, the nonequilibrium system studied in 
this paper might be accessible experimentally if the ef- 
fective temperature of the conserved magnetization com- 
ponent (s) can be maintained at a value different from 
that of the staggered magnetization which constitutes the 
nonconserved order parameter, perhaps through constant 
exposure to electromagnetic radiation. 

This work supplements earlier research that focused on 
nonequilibrium perturbations for dynamic universality 
classes which are characterized by reversible mode cou- 
plings, as relevant for second-order phase transitions in 
magnetic systems |H} (models E, G, and J, respectively 
for the critical dynamics in planar ferromagnets, isotropic 
antiferromagnets, and Heisenberg ferromagnets) as well 
as in fluids (model H for the liquid-gas transition crit- 
ical point, or more generally in binary fluids, and model E 
for the normal- to superfluid phase transition). Ref. 0] 
provides a concise summary of the results of these inves- 
tigations (including a subset of this present work) . 

We finally remark that a recent study has addressed 
a nonlocal generalization of the equilibrium relaxational 
models that allows interpolating between the scaling laws 
of models A, B, and C .20]. Intriguingly, also a 'true 
model D' scaling regime emerges in this situation, where 
both the conservation laws for the order parameter and 
the energy density remain relevant in the RG sense. 

This article is organized as follows. We start with a 
derivation of the Langevin equations of motion from the 
effective Hamiltonian for models C and D, and then pro- 
vide a brief outline of the construction of field theory 
and the perturbation expansion based on the dynamic 
action (the Janssen-De Dominicis functional |l4jp. Prior 
to describing the results of the explicit perturbation ex- 
pansion, the effect of the static nonlinear coupling of 
the order parameter to the conserved energy density is 
gauged by integrating the conserved field out of the dy- 
namic action. The implications of this, namely the re- 
duction of the isotropic nonequilibrium model D to its 
equilibrium counterpart, is discussed. Subsequently, the 
full renormalization of the vertex functions (calculated 
to one-loop order) is detailed and the expressions for the 
resulting renormalization constants are presented. From 
these Z factors we obtain the RG flow equations, and 
therefrom calculate the static and dynamic critical expo- 
nents first in equilibrium, and then successively allowing 
for isotropic violation of detailed balance in both models 
D and C. The Z factors are then adapted for the case of 
dynamical anisotropy in model C; its fixed point and crit- 
ical behavior is explained. Finally, following earlier work 
on the two-temperature model B |Tlj . the anisotropic 



nonequilibrium model D is recast into an effective two- 
temperature model D with anisotropic scaling properties. 
We conclude with a brief summary, putting our results 
into context with earlier investigations. An appendix lists 
the explicit expressions for the one-loop vertex functions. 



II. THE RELAXATIONAL MODELS C AND D 

In this section, we outline the basic model equations 
for the nonequilibrium generalization of the relaxational 
models C and D. As introduced in Ref. ^jjj], these models 
are characterized by a n-component vector order param- 
eter S = {S a },a = l,---,n, coupled to a scalar con- 
served field p. The effective Hamiltonian that describes 
their equilibrium static critical properties is the 0(n)- 
symmetric </> 4 Landau-Ginzburg- Wilson free energy in d 
space dimensions, with additional terms for the noncrit- 
ical conserved field and its coupling to the order param- 
eter. Preserving the 0{n) rotational invariance requires 
the lowest-order coupling to p to be quadratic in S. For 
models C and D, the Hamiltonian thus reads 



H[S,p}= [d d x\Y, 



^S a (x) 2 + i[VS Q (x)] 2 (2.1) 



j£^(x) 2 ^(x) 2 + fp(x)S«(x) 2 



/3=1 



Here r — (T — T°)/T° denotes the relative distance 
from the mean-field critical temperature T°, and we have 
rescaled the static energy density correlation function, 
i.e., essentially the specific heat, to unity, u and g rep- 
resent the nonlinear interaction strengths. The Hamilto- 
nian i|2.1[l determines the equilibrium probability distri- 
bution for the fields {S a } and p, 



exp(-H[S,p]/k B T) 
JV[S]V[p] eM-n[S,p]/k B T) 



(2.2) 



and furthermore provides the starting point for the field- 
theoretic static renormalization group via a series expan- 
sion in the nonlinear couplings u and g, which allows 
a systematic computation of the two independent static 
critical exponents rj and v in powers of e = 4 — d. Notice 
that the energy density fluctuations p enter the Hamilto- 
nian (|2.1|) only linearly and quadratically, and can thus 
be readily integrated out. This merely results in a shift 
of the nonlinear coupling u — > u — u — 3 g 2 . Provided the 
latter remains positive, this does not affect the RG fixed 
point u* , whence the static critical exponents are clearly 
those of the standard 0(n)-symmetric <j> model. 

We now impose Langevin dynamics on the fluctuations 
of the order parameter and the conserved field to de- 
scribe the relaxation of the system to equilibrium (at 
which the stationarity conditions 57i [S, p]/5S a = and 
6H[S, p]/Sp = hold). The purely relaxational model 
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C/D dynamics is then given by the equations of motion 

we.) = _ Mivr ^4 +r(x , i)i(2 . 3) 



dt 
dt 



Sp(x,t) 



+ r]{x,t) 



(2.4) 



Here A and D denote the relaxation coefficients of the 
order parameter and energy density, respectively (i.e., D 
is essentially the heat conductivity). The distinction be- 
tween models C and D is the value of the exponent a. For 
model C, a = corresponding to a nonconserved order 
parameter field, while for model D, a = 2, representing 
the diffusive relaxation of a conserved order parameter. 
With the effective Hamiltonian (|2.1|) . the equations of 
motion take the specific form 



ag°(x,t) 
dt 



= -A(iV) a 



(r- V 2 )S Q (x, i) + ^S Q (x,i) 



x]T^(x,i) 2 +. 9 p(x,i) S a (x,t) 



dp(x, t) 



G 



+ C(x,;), (2.5) 



dt 



p(x,i) + f£^(x,t) 5 



(2.6) 



q" . Invok- 



Upon reinstating the specific heat, note that the lin 
ear part in Eq. 1(2. 6|) implies t~ x ~ D |r| a 
ing dynamic scaling then suggests t c ~ |r 
Zp = 2 -\- a/v for the dynamic exponent of the energy 
density. We shall see that this relation in fact holds only 
for n < 4, or more precisely, when a > 0, for otherwise 
the two Langevin equations decouple at criticality. 

In the above stochastic equations of motion 1)2.5(1 and 
(|2.ti|) . ( a and 77 represent the stochastic forces ('noise') for 
the order parameter and the conserved field, respectively. 
We assume a Gaussian distribution for these fast vari- 
ables with a vanishing temporal average, (£ a ) = = (77). 
Their second moments then take the functional form 

(C Q (x,i)C ,3 (x',t')> = 2A (iV) a <5(x - x') 6(t - t') 8 ap , 

(2.7) 

(77(x,i)r7(x',0) = -2D\7 2 5{x-x')5{t-t'). (2-8) 

In thermal equilibrium, Einstein relations connect 
the relaxation coefficients with the corresponding noise 
strengths according to A = A k-gT and D = D fceT, 
where T is the temperature of the heat bath in contact 
with the system. The detailed balance conditions im- 
plicit in these relations ensure that the system relaxes to 
the equilibrium probability distribution 12.21) in the limit 
t — ► 00. More generally, we can identify A/ A = k^Tg and 
D/D = kftTp as the temperatures of the heat baths cou- 
pling to the order parameter and the conserved density, 
respectively, with their ratio given by 



e 



In 



D A 
D A 



(2.9) 



This new degree of freedom 8 describes the extent to 
which the equilibrium condition is violated; detailed bal- 
ance clearly holds for 8 = 1 , while for 8 < 1 energy flows 
from the order parameter heat bath to the conserved den- 
sity heat bath, and vice versa for 8 > 1. Since we are 
interested in the behavior near the critical point Ts ~ T c , 
8 essentially measures the temperature of the conserved 
density heat bath T p , in units of the critical tempera- 
ture T c . In the critical regime, we will be interested in 
calculating the dynamic exponents zg and z p that de- 
scribes the critical slowing down for the order parameter 
and the conserved energy density. The one-loop RG flow 
equations for this case of isotropic detailed balance vi- 
olation are derived in Sec. IIIII The assumed functional 
form for the noise correlations ((2.7(1 , ((2.8(1 also enables us 
to impose a spatially anisotropic form of detailed balance 
violation, as described in Sec. IIVI 

We close this general introduction with a brief outline 
of how a field theory representation can be constructed 
from general Langevin-type equations of the form 



d#*(x,t) 
dt 



ini> a }](x,t) + C Q (x,i) . (2-10) 



with (£") = 0, and the general noise correlations 
(C"(x,t)C /3 (x',i')) =2L a 5(x-x.')5(t-t')5 af} . (2.11) 

This form of the white noise may be inferred from a Gaus- 
sian distribution for the stochastic forces 



IU[{C a }] oc exp 



-i Jd d x JdtJ2C(L a )-\ a 



(2.12) 

Eliminating £ Q via Eq. ((2.1(J(I then immediately yields 
the desired probability distribution for the fields i[i a , 

W{{( a }}V[{( a }} = P[{^ a }]V[{i) a }] oc e G[{ ^ }] V[{^ a }] , 

(2.13) 

with the Gnsager-Machlup functional 



cur}] 



d d x I dt 



x(L c 



^\ dt 

diji a 



K a [{r}\ 

-K a [{r}]) ■ (2.14) 



From this functional, one can already construct a per- 
turbation expansion for the correlation functions of the 
fields ip a ; however, since the inverse of the Onsager co- 
efficient L a is singular for the conserved quantities, and 
furthermore high nonlinearities oc K a [{ijj a }}' 2 appear, it 
is convenient to introduce Martin-Siggia-Rose auxiliary 
fields via a Gaussian transformation to partially linearize 
the above functional. This leads to 



par}] cx j mr}] exp(-^[w Q }, {r}]) (2.15) 

with the Janssen-De Dominicis functional 0] 



A[{r},{r}} 



J d d x JdtJ2 
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K a [{r}} 



(2.16) 



Eq. I|2.16(l will provide the starting point for our discus- 
sion of the nonequilibrium dynamics of models C/D. In 
Sec. IIIII we will use the corresponding Janssen-De Do- 
minicis functional for the construction of dynamic per- 
turbation theory, and therefrom infer the one-loop RG 
flow equations, first in equilibrium and then with broken 
detailed balance. Subsequently in Sec. lIVI we will repeat 
the procedure for the models with anisotropic detailed 
balance violation. 



III. THE ISOTROPIC NONEQUILIBRIUM 
MODELS C AND D 



In this section, we will study the nonequilibrium crit- 
ical properties of the relaxational models C and D with 
isotropic violation of detailed balance. Along the way, 
we shall also recover the equilibrium critical exponents. 
The field theory is constructed as outlined in the pre- 
ceding Sec. [HI and a perturbation series expansion in 
the relevant nonlinear couplings cx u and g 2 is developed 
for the one-particle irreducible vertex functions, explic- 
itly here to one-loop order. The subsequent renormaliza- 
tion constitutes a straightforward generalization of the 
equilibrium renormalization scheme, see Ref. |T6j . From 
the renormalization constants (Z factors) that render 
the field theory finite in the ultraviolet (UV), we derive 
the RG flow functions which enter the Callan-Symanzik 
equation. This partial differential equation describes the 
behavior of the correlation functions under scale trans- 
formations. In the vicinity of an RG fixed point, the 
theory becomes scale-invariant and the information from 
the UV behavior can be employed to access the physi- 
cally interesting power laws governing the infrared (IR) 
regime at the critical point (r cx T — T c — > 0) , for long 
wavelengths (q — > 0) and at low frequencies (cj — ► 0). 



A. Dynamic field theory for models C and D 

As a first step, we translate the Langevin equations 
(|2.5() and l|2.6(l , with the noise correlations l|2.7[l and 12.811 , 
to a dynamic field theory 0, . This results in a prob- 
ability distribution for the dynamic fields S and p 

P[S,p]« Jv[{iS a }} fv\ip\ exp(-A[a,S,p,p], (3.1) 

with the statistical weight given by the Janssen-De Do- 
minicis functional A = Ah&r +-4rel +-<4cd- The harmonic 
part, in terms of the original dynamic fields S a and p, 
and the corresponding auxiliary fields S a and p, reads 



Aar[S,S,p,p] = Jd d xjdt (^S Q (x,i) 



dS a (x,t) 

at 



-A (iV) Q (r - V 2 ) 5 Q (x, t) - A (iV) a S* a (x, t) 
~dp(x,t) 



(3.2) 



+p(x, t) 



Of 



DV' p(x,t) + DVp(x, t) 



while the static nonlincarity leads to a relaxation vertex, 

Ael[S, S] = ^Jd d X Jdt £ S Q (X, t) 



a,/3 



A(zV) a ^(x,i) 2 



S a (x,t) , (3.3) 



and the coupling between the order parameter and the 
conserved density generates the model C/D vertices 



Ad[S, S, p, p]=g J d d x JdtJ2 



xp(x, t) S Q (x, t) - p(x, t) D V 2 - S a (x, tf 



S a (x,t)X(iV) a 
(3.4) 



Before we proceed to develop the perturbation expan- 
sion based on the above dynamic functional, we can try to 
gauge the relevance of the nonequilibrium parameter 0, 
as defined in Eq. (|2.9I) . by integrating out the conserved 
density p from the action. Denoting those terms in the 
total dynamic action A[S,S,p,p] that involve only the 
order parameter and the corresponding auxiliary fields 
as A[S, S], and subtracting this part, we are, in Fourier 
space, left with 



A[p,p] = A[S,S,p,p]-A[S,S] 
d d q fduj f 



(3.5) 

H^ + £>q 2 )p(q,w) 



,Dq 2 p(q, W ) + Z?q 2 |5 2 ( q , W ) 



{2n) d J 2tt 



+p(q,u;)\q a g[S-S](-q, -lo) 
where we have introduced the composite operators 

s 2 (q,")= (3-6) 

v ' a 

[5.S](q,w)= (3.7) 
d d p f dv 



bJ — V) 



as Fourier convolutions. 

The path integral over the fields p and p now takes the 
form, in matrix notation, 

V\ip\Jv[p] exp(-A[p,p]) (3.8) 

= Y[Jv[ix 1 )V[x 2 ]exp^-^x T Ax~b T x^ , 
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with the vectors 



p(q.w) 


, b = 


. p(q^) _ 





Dq 2 ffS 2 (q,w)/2 
Xq a g[S-S](q^) 



and the Hermitean matrix 
A = 



-2Dq 2 -iuj + Dq 2 
iw + D q 2 



(3.9) 



(3.10) 



After the linear transformation y = x+A 1 6, the integral 
becomes 



\y T Ay\ exp [\b T A-H 



V[iyi\D{y 2 ] exp 

where the entries of the inverse matrix 

A- 1 - 



(3.11) 



Jiuj + Dcf)- 1 

(-iuj + Dq 2 )- 1 2Dq 2 /{uj 2 + L> 2 q 4 ) 



(3.12) 

actually represent the propagators for the scalar con- 
served held p. Upon performing the Gaussian integration 
over the y fields, the p and p fields are integrated out to 
yield the effective dynamic functional 



■Aeff[S,S] = .A[S,S 
d d q 



(3.13) 



(2tt) 



^ AqV [S-SK-q.-u;) 

Z7T 



S 2 (q,uj)/2 X Dq a [S ■ S}(q,u)/Dq 2 



l-iu/Dq 2 D \ + {u)/Dq 2 f 
We now define the parameter 

w = D/X , 



(3.14) 



which essentially measures the ratio of relaxation times 
of the order parameter and the conserved field, i.e., 
w ~ Tg/T p . For model D (a = 2) at criticality, the relax- 
ation time of the conserved order parameter field is much 
longer (since dS a /dt ~ q 4 ) compared to that of the con- 
served field (since dp/dt ~ q 2 ), so that w — > oo as q — > 0. 
Hence the second term in the brackets in the^ above ef- 
fective functional (which contains the ratio D/D) van- 
ishes asymptotically. Consequently, the noncquilibrium 
parameter 9 disappears from the effective field theory 
(|3.13|) entirely, and a simple rescaling of the nonlinear 
couplings u and g reduces model D with isotropic de- 
tailed balance violation to its equilibrium counterpart. 
This remarkable result will be borne out in the explicit 
one-loop perturbation theory as well, see Sec. IIIIDI 



B. Perturbation theory and renormalization 

1. Elements of dynamic perturbation theory 



C and D. The harmonic part (|3.2() defines the (bare) 
propagators of the field theory, while the perturbation 
expansion is performed in terms of the nonlinear vertices 
(|3.3|) and l|3.4|l . Note that the existence of the irreversible 
forces l|3.4fl does not show up in dynamic mean-field the- 
ory (van Hove theory) at all, which is based on the har- 
monic action <|3.2|l only. 

We can now construct the perturbation expansion for 
all possible correlation functions of the dynamic and aux- 
iliary fields, to be computed with the statistical weight 
exp(— A[S, S, p, p\), as well as for the associated vertex 
functions given by the one-particle irreducible Feynman 
diagrams. A straightforward scaling analysis yields that 
the upper critical dimension of this model is d c = 4 
for the relaxational vertices l|3.3|) and l|3.4|) . Therefore, 
for d < 4 the perturbation theory will be IR-singular, 
and nontrivial critical exponents will result, whereas for 
d > 4 the perturbation theory contains UV divergences. 
In order to renormalize the field theory in the ultravio- 
let, it suffices to render all the nonvanishing two-, three-, 
and four-point functions finite by introducing multiplica- 
tive renormalization constants (in addition to an additive 
renormalization that amounts to a fluctuation-induced 
shift of the critical temperature). This is achieved by 
demanding the renormalizcd vertex functions, or appro- 
priate momentum and frequency derivatives thereof, to 
be finite when the fluctuation integrals are taken at a con- 
veniently chosen normalization point p, well outside the 
IR regime. Notice that p defines an intrinsic momentum 
scale of the renormalized theory. The Callan-Symanzik 
equations can subsequently be used to explore the de- 
pendence of the renormalizcd vertex functions on p, and 
thereby obtain information on the scaling behavior of the 
dynamic correlation and response functions. 

The Gaussian (zeroth-order) propagators 

G| Q5 ,(q,-) = r|^(-q,-.)-\ (3.15) 
G&(q, W ) = ry-q,-*)- 1 (3.16) 

and vertices which are the starting point for perturbation 
theory are 



T~ „ (q, w 

r9~(q,w 
pp 



^ + Aq a {r + q 2 )]S al3 . 
iui + Dq 2 , 

2Xq a S af3 , 
2Dq 2 , 
-Xq a gS^, 

-\Dq 2 g5^ , 

\ a U 

-Aq°- . 
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(3.17) 
(3.18) 
(3.19) 
(3.20) 
(3.21) 

(3.22) 
(3.23) 



We first detail the dynamic field theory for the case 
of isotropic detailed balance violation for both models 



In addition to the previously introduced ratio of relax- 
ation times w and the nonequilibrium parameter 0, we 
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define for convenience the rescaled static couplings 



u = — u 

X 



~2 ^2 

A 



(3.24) 



Recall that is a measure of the extent to which detailed 
balance is violated, since for = 1 a straightforward 
rescaling of the couplings reduces the dynamic functional 
(|3.2(l - (|3.4|l to the equilibrium case. 



2. Vertex function renormalization 

The explicit expressions for the relevant vertex func- 
tions to one-loop order in the perturbation expansion 
are given in appendix ^ The ultraviolet-divergent 
derivatives of the two-, three-, and four-point vertex 
functions that require multiplicative renormalization are 
<Vr~ s (q,0) | g= o, 5 9 2 +<I r~ s (q,0) | g= o, d^T~ s (0,uj) | w=0 , 



<9 g 2r~ p (q,o) \ q=0 , a w r~ p (o,tj) | w = , 9 g «r— (q,o) | 9=0 , 

<9 9 2r~(q,0) | g= o, d qa r~ Sn (q,0) \ q=0 , d q 2T~ ss (q,0) \ q = , 



SSp 



and d q aT~ sss (q, 0) | g= o- The quadratic divergence in 
the first of these will be taken care of by the T c shift r c . 
The remainder as well as all the other expressions are 
logarithmically divergent at the upper critical dimension 
d c = 4. We thus require ten multiplicative renormal- 
izations in all, which we take to define the renormalized 
counterparts of the fields S*, 5, p, and p, and of the pa- 
rameters -D, D, A, A, it, g, and r = r— r c , which represents 
the temperature distance from the true critical tempera- 
ture. The renormalized quantities are defined through 



5a ^1/2 Qa qc 



ryl/2 ~ 
PR = Z p PR 



ryl/2 ~ 

ZJ p. 



(3.25) 
(3.26) 

t r = Z t t p~ z , X R = Z x A , X R = Z~\ , (3.27) 
D R = Z D D , D R = Z~ D , (3.28) 
ur = Z u u A d pT e , g% = Z g g 2 A d p~ e , (3.29) 



where in standard notation 



and A n 



r(3- d/2) 

2<J-l 7r d/2 



(3.30) 



The loop integrals are evaluated in the dimensional regu- 
larization scheme, and we choose the renormalized mass 
tr = 1 as our normalization point (i.e., r = p?Z~ x w p 2 , 
to lowest order). Notice that as we have thus defined 
eleven renormalization constants (Z factors), but actu- 
ally only need ten multiplicative renormalizations, we 
shall have the freedom to choose, say, Xr — Xr. As 
A = A can be achieved through simple rescaling in the 
unrenormalized theory, this implies the choice Z~ = Z\. 
In addition, the structure of the perturbation series im- 
plies nontrivial additional identities between the renor- 
malization constants, as we shall see below. 



We now proceed to compute the renormalization fac- 
tors by absorbing the UV divergences of the loop inte- 
grals following the minimal subtraction prescription. All 
subsequent explicit one-loop results are for the case of 
model C (a = 0). However, as we have seen at the end 
of Sec. 1111 Al the separation of relaxation time scales for 
the order parameter and the conserved density in model 
D leads to w — > 00 in the asymptotic limit. It turns out 
that taking this limit for the model C results precisely 
yields the Z factors for model D. An independent calcu- 
lation of the Z factors from the model D vertex functions 
confirms the validity of this simple limit procedure. 

First, we employ the criticality condition x(0,0) _1 = 
r~ s (0,0)/A = at the true critical point r — r c , and 
solve for the fluctuation- induced T r shift, 



Ta = -- 



(u-Sg 2 ) 



-9 2 



1-6 



l + io J p r c /(l + w) +p 2 



(3.31) 



We may then reparametrize r~ s (0,0) in terms of r 
r — r c , which amounts to an additive renormalization, 



r~„(0,0) = At 



1 - 



-g 2 



l-e 



(u-Sg 2 ) 
1 



1 



(1 + w) 2 J p p 2 [r/(l + w) + p 



P P 2 (t + P 2 ) 
(3.32) 



Writing this result in terms of renormalized quantities, 
and evaluating the integrals at the normalization point 
t = p 2 in dimensional regularization, we obtain the fol- 
lowing expression for the product of the Z factors 



{Z~Z S ) 1/2 Z x Z r = 



(3.33) 



G 



(u-3g 2 )+g 2 



i-e 



(l + w) 2 



A d p~ 



Next, expanding the integrands in the expression for 
r~ s ,(q, oj) to order q 2 to obtain the renormalization fac- 
tor for the relaxation rate A, we find 



d_ 
dq 



+ ■ 



(3.34) 



g 2 1 



4 1 
f (1- 



e 



6)(1 



(1 + w) 3 7 p [r/(l + w)+p 



2i:s 



whence 



{Z~Z S ) 1/2 Z X = 1- 



g 2 i-e 



4 1 



1 - 



(1-w) 



2 1 



(1 



A d p~ 



Another product of Z factors is obtained from 




d(ioj) T s> 



j=0 



(3.35) 



(3.36) 
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= 1 



which gives 

{Z~Z S ) 1/2 = 1 



1 + w J p (r + p 2 ) [r/(l + w) + p 2 

i-e y i 

(1 + W ) 2 ^^/(l + ^+p 2 ] 2 : 



g 2 

1 + 10 



1-9 



(3.37) 



Now consider the two-point vertex function T~^(q, w) 
for the conserved density. Because any loop diagram 
for this quantity necessarily involves a pSS vertex for 
its outgoing leg, we see that taking q — > results in 
T~ p (0,Lu) = ilu to all orders in the perturbation expan- 
sion. As a consequence, 



Z~Z n 



1 . 



Upon absorbing the logarithmic divergence of 

1 



r~( q ,o) = iv 



(r + p 



2 V2 



into Zjj, we arrive at 

Z D = 1 



7:2 



(3.38) 



(3.39) 



(3.40) 



The vertex function T~(q, ui) is actually UV-finite to all 
orders. Again from the momentum dependence of the 



pSS vertex, d q ? Y — (q, 0) | g =o= —2D, whence 

Z~Z~ ee 1 , (3.41) 

and with Eq. (|3.38|l therefore 

Z~ ee Z p . (3.42) 

The remaining logarithmically divergent two-point func- 
tion 



r^(o,o) 



(3.43) 



-2 A 



e 



1 + w J p (t + p 2 ) [t/(1 + w) + p 2 ] 



yields the relation 



e A dl r 



1 + w 



(3.44) 



At last, from the rather lengthy one-loop results (|A5|I - 
(|A7|) for the three- and four-point vertex functions we 
deduce, respectively, 



ZsiZ-Zg) 1 ' 2 Z D = 

1 n + 2 „ A d fi- 
1 — u 



~2 

9 



1-6 



1 + u> 



{ZgZs Z p Zg) 1 / 2 Z\ — 



(3.45) 



(3.46) 



6 

a/2 



2 ~ A di i 
■ u 



+ ■ 



i-e 



(l+w) 2 



A d p 



{Z-ZsY 1 ' Z S Z X Z U = \- 



A 



(in 



(i 



+60" 



1 - 

u 



(l-9)(2 + w) 
2(1 + w) 2 



e 

A d p 



1 - 



i -e 



l + w 



A d fi 



(3.47) 



Upon factoring from the above Z factor products 
(|3~33|) . E2S|), (|3~4T)|) . and (JOSHSSZIj the fol- 

lowing results are obtained: 



n + 2 
Z r = 1 — u 



Z\ = l- 



6 e 
n + 2 1-9 



l + w 

2 



1 - 



(1+W) 3 

1-9 



.4 



(1 + w) 2 



e 

AdV 



, (3.48) 
, (3.49) 
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n _ , _ n] w(2 + w) 
2 1 1 {l + wf 



Ad\±~ 



Z s (Z- 1 Z g ) 1 / 2 = l- 



2 _ A d \i~ 
— u 



+9 1 



1-9 



l + w 



A d fJ,~ 



(z P z g y/ 2 = i 



2 _ A d p 
— u 



-a 2 



e 

1-9 



(l + w) 3 



A d p,~ 



Z s Z a = l- 

+9 
Zs Z u = 1 
6 



2 _ Adp 
■ u 



3 

4 _ _ 2 + 2w + w 2 



(3.50) 
(3.51) 

(3.52) 
(3.53) 



2 (l + w) 3 

n + 8 „ Ad/j, 



A d p 



+9* 





6 ° 


e 




1 


- 9 /3(2 


+ w) 




1 


+ w V 


1 - 


f w 








1 - 


1 - 


9 




u 


14 


w 



w 



(l + w) 2 
A d [i-' L 



A&\x~ 



(3.54) 



supp lemen ting Eqs. ESHl, ^3QJ, (EIJ), 

and JSIllJl' 

For model D (a = 2), the ratio of relaxation times w 
constitutes a relevant parameter in the RG sense, whence 
w — * oo asymptotically. Obviously, this leads to marked 
simplifications in the above expressions for the renormal- 
ization constants. In addition, as a consequence of the or- 
der parameter conservation law and the ensuing momen- 
tum dependence of the vertices, we have r~ s (0,w) = iuj 

and <9 g 2 r~~(q, 0) | 9 =q= —2 A to all orders in pertur- 
bation theory, which implies the relations Z~Zg = 1, 
Z~ Z-r = l, and thus Zr = Zs, which also follow to one- 
loop order from Eqs. i|3.37[) and l|3.44|) . respectively. 
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Callan-Symanzik and RG flow equations 



C. Model C and D equilibrium critical exponents 



By means of the above renormalization constants, we 
can now write down the Callan-Symanzik RG equations 
for the vertex functions and the dynamic susceptibilities, 
which describe the dependence on the renormalization 
scale fj,, and thus on the renormalized couplings. These 
RG equations connect the asymptotic theory, where the 
IR singularities become manifest, with a region in param- 
eter space where the loop integrals are finite and ordinary 
'naive' perturbation expansion is applicable. They follow 
from the observation that the 'bare' vertex functions do 
not depend on the renormalization scale /j,, 



dfi 



r 



S m p"5 r p s 



({q, W };{a}) = 



(3.55) 



where {a} represents the parameter set u, g 2 , D, D, A, 
A, and r. Replacing the bare parameters and fields in 
Eq. (|3.55() with the renormalized ones, we find the fol- 
lowing partial differential equations for the renormalized 
vertex functions: 



d ^ d 
V + 1^ a R - — 
on f-C oa R 



777 

-77 11 



— 7~- 
2 'p 



2 13 



R 



({q,u;};{a fl }) = 0. (3.56) 



Here, we have introduced Wilson's flow functions 



d_ 



for the fields and 



lnZ~ 



lnZ~ 



is 



d_ 



lnZ„ 



la = M^- 
0(1 



In 25 



(3.57) 
(3.58) 

(3.59) 



for the different parameters (the subscript '0' indicates 
that the renormalized fields and parameters are to be ex- 
pressed in terms of their bare counterparts prior to taking 
the derivatives with respect to the momentum scale fi). 

The Callan-Symanzik equations l|3.56[l are solved via 
the method of characteristics, introducing fi(£) = fj£, 
where £ is a real continuous parameter. This defines run- 
ning couplings as the solutions to the first-order differen- 
tial RG flow equations 

da{£) 



I- 



d£ 



l a (£)a(£) , with a(l)=a R 



(3.60) 



The solutions of the partial differential equations (|3.56() 
then read 



(Mq,^};{ai?J) = 



(3.61) 



exp 



e d£' r 



■7s(0 



■lp( 



(p£,{q,u};{a(£)}) 



To begin the analysis of the RG flow equations, we re- 
cover the critical exponents for the equilibrium models 
C and D (see the original Ref. and Ref. 0] for the 
corresponding field the ory; the two-loop analysis was re- 
cently clarified in Ref. (17|). Upon removing the effects 
of the nonequilibrium perturbation by setting A = A and 
D = D (we may put k^T « k^T c = 1), whence 8 = 1, 
7t = 7_t, and g 2 = g 2 in the previous expressions for the Z 
factors, we obtain the renormalization constants for the 
equilibrium case to one-loop order: 



Z p =Z~ 
p 

Zs = 1 ■ 



Zr 



n g 2 AdfT 
2 e 



1 g 2 Adfi 



(3.62) 
(3.63) 

(3.64) 



zz 1 ' 2 = Z X = 1- 
s i + w e 

77 + 2 uA di i- € 77 + 2 g 2 A d ^ 
Zr = l~— 6 — + — —,(3.65) 

n + 2 uA d n- e n + 4 g 2 A df i~ e 
1 1 , (3.66) 



Z„ 



3 e 

77 + 8 uA d fL~ 



6 



6 1 



e 

g 2 \ g 2 A d n~ 



(3.67) 



In this equilibrium system, there exists a fluctuation- 
dissipation relation that relates the imaginary part of 
the dynamic order parameter susceptibility x(q, ui) to 
the Fourier transform of the dynamic correlation func- 
tion C{x-x',t-t')6 a P = (S a U,t) SP{x.',t')): 

C(q,oj) = ^- Im X (q,w) . (3.68) 

LO 

These quantities are connected to the two-point ver- 
tex functions via C(q,w) = -r~~(q, w)/|r~ g (q, u>)\ 2 , 
and x(q, w) = Aq a /T gg (— q, — u)\ thus the fluctuation- 
dissipation theorem results in the following relation be- 
tween the two-point vertex functions 



2A 

r~~(q,tj) = — Imr~ (q,w) 



(3.69) 



The same identity must hold in the renormalized theory. 
Consequently, 



Zx = 



(Z S /Z~) 



1/2 



and in the same manner for the conserved field, 

1/2 



J D 



(z P /z~) 



(3.70) 



(3.71) 



Both relations are indeed fulfilled by the above explicit 
one-loop results. Through taking logarithmic derivatives 
with respect to the normalization scale /x, one finds that 
the equilibrium fluctuation-dissipation theorem implies 



2 7a = IS - 7o and 2 j D = 7 P - 7~ 



(3.72) 
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The explicit RG flow functions derived from the one- 
loop renormalization constants become 



7s 

7p : 



„, 7J = _| = ^L_ 

2 1 + wr 



"7~ = ID 



9r 



7r = -2 - 



6 



UR 



(3.73) 
(3.74) 
(3.75) 



In the result for the single nontrivial static RG flow func- 
tion i|3.75|) to one-loop order, ur — ur — 3 g R represents 
the shifted coupling that also results from directly inte- 
grating out the scalar density p, see Sec.|n] Indeed, from 
Eqs. (|3.66[> and l|3.67[l we infer its Z factor 



n + 8 uAju 



6 



(3.76) 



which along with Eq. I|3.75|l is just the standard one-loop 
result for the 0(n)-symmetric 4 theory. Note that for 
model D, taking w R — > oo in Eq. I|3.73|) yields 7^ = 0. 

We are now in a position to study the scaling behavior 
in the vicinity of the various RG fixed points which are 
given by the zeros of the RG functions 



0a = ladR = A* 



d_ 

Of! 



CLR 



(3.77) 



for the nonlinear couplings u and g 2 as well as the re- 
laxation rate ratio w, W = wr (70 — 7a). By means of 
Eqs. (|3~52|) . £HS|), (|3U6) . and JTZfJ, we find 



0u 
09 
0w 



Ur 



9r 



n + 8 _ 
6 

n + 2 



wr 9r 



— UR + -9R 



1 



1 + WR, 



(3.78) 
(3.79) 
(3.80) 



For model C, the flow function W yields three fixed 
points, provided g* 2 > 0, namely Wq = 0, = (2/n) — 1 
(which is positive for < n < 2), and w* D — 00. Stability 
requires that 



9l 



1 



(1 + W R f 



(3.81) 



be positive at the fixed point. Consequently, for n = 1 
we find that = 1 is stable, whereas Wq = for n > 2. 
Recall that w^, = 00 corresponds to model D; this fixed 
point is unstable in model C for all values of n. 

For the static coupling ur we find the following zeros of 
0u : the Gaussian fixed point Uq = 0, and the Heisenberg 
fixed point u* H — 6 e/(n + 8). Inserting these in turn into 
the flow function [3 g , we obtain: <?q 2 = and g\ 2 — 2 e/n 
corresponding to Uq = 0; g^ 2 = and g* c 2 = 2 (4 — 
n) e/n (n+8) corresponding to u* H . The stability of these 
four fixed points in the {ur^r) plane depends on the 



spatial dimension d and the number of order parameter 
components n. Checking for positivity of the eigenvalues 
of the stability matrix with the entries dfiu/duR = — e + 
(n + 8) ur/3, dPuldg 2 R = 0, d(3 g /du R = (n + 2) g 2 R /3, 
and d0 g /dg 2 R = — e + (n + 2) u R /3 + ng R at these RG 
function zeros, we find that for e > (d < d c = 4) the 
only stable fixed points to one-loop order are: [u^,g^ 2 ], 
stable for < n < 4, and [u* Hl gl 2 ] 7 stable for n > 4. 

The solutions of the RG equations (|3.61() yield for the 
order parameter susceptibility and correlation function in 
the vicinity of an IR-stable fixed point the scaling laws 



2+a+ 7 * 



) , (3.82) 



C(r, q,w) 



-a-7x-7s 



Setting to = 
exponents 



in Eq. 



s C{Tq^,Uj/q 2+a+ ^) . 

(3.83) 

we identify the static critical 



V = -7s = , 
- 1 = -7* = 2 - 



(3.84) 
(3.85) 



with their one-loop values computed at the Heisenberg 
fixed point u* H . The standard hyperscaling relation then 
gives for the critical exponent of the specific heat 



a 



2-dv 



Therefore, we may rewrite 



9c 



2 (n + 8) 

2 a 
y = — — 

n v 



(3.86) 



(3.87) 



to this order in e = 4 — d. The dynamical critical ex- 
ponents zs and z p that describe the divergence of the 
characteristic relaxation times for the order parameter 
and the conserved density, respectively, are given by 



zs 



a + 7A 
1*D ■ 



(3.88) 
(3.89) 



For model C (a = 0) with nonconserved order pa- 
rameter, we thus have three equilibrium scaling regimes 
[T5I fit!, lri| : (a) In the first regime with n — 1 (Ising 
symmetry), the stable critical one-loop fixed point is 



•h 



2c 



9c 



2e 

Y 



1 



(3.90) 



This describes a strong scaling regime where the dynamic 
exponents for the order parameter and the conserved field 
are identical, 



z s = z p = 2+ - = 2+ - 
3 v 



(3.91) 



(b) In the second regime with 2 < n < 4, the stable 
critical fixed point becomes 



' H 



6e 
n + 8 



9c 



2 (4 - n) e 
n(n + 8) 







(3.92) 
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leading to weak dynamic scaling with 

2 (4 - n) e „ 2 a „ a . . 

z s = 2 + 7 , g , = 2 + -- <z p = 2 + -. 3.93 
n (n + 8) rt f f 



(c) Lastly, for n > 4, 
6e 



5o = , Wo = 



(3.94) 



and consequently z$ = z p — 2 take on their mean-field 
values to one-loop order. More generally, for a < the 
order parameter and conserved energy density dynamics 
decouple at criticality, which implies purely model A dy- 
namics for the order parameter, and uncritical diffusive 
relaxation for the conserved mode, i.e., 

z s = 2 + cr], z p = 2, (3.95) 

withc = 6 1nf-l + 0(e). 

For model D, the conserved order parameter (a = 2) 
always relaxes much slower than the also conserved, but 
noncritical energy density near the phase transition, and 
consequently w — > oo. Notice that the identity Z\ = Zs 
implies 7a = 7s, and hence the model B scaling relation 
zs = 4 — 77 holds. We now have only two different regimes, 
with the conserved field either influenced by the critical 
variable, or not. For n < 4 {a > 0) 



•H 



whence 



6e 
n + 8 



9c 



2 _ 2(4-n)e 
n (n + 8) 



, w D = 00 



zs = 4 - 77 , z p = 2 + - , 

v 



whereas for n > 4 (a < 0) 
6e ,9 



9o 







w D = 00 , 



(3.96) 



(3.97) 



(3.98) 



with the decoupled model B dynamics described by 

Z S EE 4 - 7] , Zp EE 2 . (3.99) 

D. Isotropic detailed balance violation in model D 

We start with a particularly simple case of our various 
nonequilibrium systems, namely that of model D sub- 
ject to isotropic detailed balance violation. Previously, 
we saw upon integrating out the conserved field from the 
dynamic action and taking the asymptotic limit of the 
ratio of relaxation times, i.e., w = D/X — * 00, that the 
nonequilibrium parameter G drops out of the field theory 
entirely This is also explicitly seen in the perturbation 
expansions for the vertex functions. For example, con- 
sider the two point function 



r~ s (q,0) = Aq : 
+9 2 (l-&) 



r + 



(u-3g 2 ) 
(§+P) 2 



1 



(f +P) 2 [r + (f + P ) 2 ]+w(f-p)2 



c.f. Eq. (HJ in appendix 1X1 We see that the second term 
which involves vanishes in the limit w — > 00. This 
is generally true of all the other vertex functions also. 
Therefore, all the effects of the nonequilibrium perturba- 
tion are eliminated. Consequently, the renormalization 
factors then become identical to those of the pure equi- 
librium model (|3.fi2|l - H3.67J) . with however, the rescaled 
nonlinear couplings u u and g 2 — > g 2 . Therefore, the 
critical fixed points and exponents are also accordingly 
identical to those of the equilibrium model D. This leads 
us to the distinct statement that a conserved order pa- 
rameter subject to an isotropic nonequilibrium perturba- 
tion does not display any novel dynamic critical behav- 
ior even when quadratically coupled to a conserved scalar 
density. We shall, however, see in Sec. iTVCl that when 
model D is subject to an anisotropic 'dynamical' noise, 
drastic effects may emerge in this system. 



E. Isotropic detailed balance violation in model C 

1. Renormalization and one-loop RG flow functions 

For model C (a = 0), w < 00 at the stable equi- 
librium fixed points, so the nonequilibrium parameter 
does not disappear from the asymptotic theory. As 
mentioned before, a simple rescaling of the fields and 
coupling constants allows setting the relaxation rate and 
noise strength of the order parameter equal, A = A (with 
k B T s = 1 here), whence = D/D. With Z~ = Z X , the 
ratio of Eqs. (|3.44|) and (|3.49|l gives Z~, and subsequently 
by means of Eqs. (|3~37|) . (|3~5U|) - (|3~P|l we arrive at 



Z s 

Z n : 



f) 2 



n 2 (1 - 0) w 
2 ~ {1 + w) 2 

2 (1 - g) ™ 2 A d fi~ 
(1 + w) 3 e 
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(3.100) 
(3.101) 
(3.102) 
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2 „ A d \i 
— u 
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n + 4 2(1-0) 
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n + 8 _ A d \jT 
Z u = \ — u ■ 



6 e 

2(1-0) 
l + w 



0- 




(i 
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u 




1 




l + w) 





w 



1-0 



l + w 

A d pL- c 



(3.103) 

A d ti~ e 
e 

A d (j.~ e 



(3.104) 



From those renormalization constants, we infer the RG 
flow functions 



7S 



_a (1-Qr)w 2 r 



'9b. 



(1 + w r ) 3 



(3.105) 
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7c = 



9% 



1 + Wr 



9% 



1 + R - 

1 - 



(i + ^fl) 2 
i-e R 



7p = 

7r : 



-7r = 7s = fl| 



(1 + ^) 2 
n 2(1-9^)^ 
2 



~2 



(3.106) 

(3.107) 

(3.108) 
(3.109) 



n + 2 



U R ~ 9r 



n + 2 



I-Qr ' 
(1 + Wfl) 3 . 

(3.110) 



and the four coupled RG (3 functions 

i-e R 



i 



-2& 



1 + Wj{ 

1_6r 

1 + WJ.R 

n + 2 



n + 8 _ 
+ w/i -e H — u R 



6 



1 



1 + W_R 



, (3.111) 
(3.112) 



n + 4 1 - 6^ 
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1 + W7f V (1 + wt?J 2 
n 1 1 - O fl 

2 ~ l + w R {1+wr) 3 



, (3.113) 



/3e = e J? (7~-7z 5 ) = -2g 2 , 



WrQr(1-Qr) 



(3.114) 



2. RG fixed points and their stability 



For wr g% > 0, Eq. (|3.114|l yields three RG fixed points 
for the nonequilibrium parameter Q, namely 0g = 0, 
oo. But from 



6* q = 1, and e 



dfo ^w R {l-2Q R ) 
— — z 9r ■ 



(3.115) 



we infer that only the equilibrium fixed point 0* q = 1 is 
stable. This implies that detailed balance is effectively 
restored at the phase transition in this situation, and the 
asymptotic critical behavior is that of the equilibrium 
model C with the exponents given in Sec. lIII 01 It is how- 
ever instructive to investigate the possible existence of 
genuine nonequilibrium fixed points, as they might influ- 
ence the scaling behavior in transient crossover regimes. 

We begin with 6q = 0. The RG function for the 
time scale ratio wr at Or — reads 



A, 



W R 9r 



with the derivative 



d(3 w 
dwR 



9R 



n _ Wr (2 + Wr) 

2 {l + w R f 



n Wr (4 + Wr) 

2 



(1 + wr) 4 



(3.116) 



(3.117) 



Since the maximum value of the second term in the brack- 
ets of Eq. I|3. 11611 at w R = a/3 - 1 is w 0.385 < n/2 for 
any n > 1, the only fixed points are Wq = and w* D = oo, 
of which the former is stable. At 8q = = Wq, we find 



7s = 7- 
7 P = -7- 

1* = -2 • 



7a =1-1 



Id = lj 



6 



2 9 



2 9 



with u* and g* 2 denoting the zeros of 



% = 9r 



n + 8 
6 

n + 2 



UR 



(3.118) 
(3.119) 

(3.120) 

(3.121) 
(3.122) 



Note the striking similarity with the equilibrium (3 func- 
tions H3.78fl and l|3.79|) . yet with the crucial sign change 
in Eq. (|3.122|) . and the fact that the anomalous dimen- 
sions at the fixed point satisfy the relations l|3.72|l that 
would be imposed by a fluctuation-dissipation theorem. 
Upon inserting the stable (for e > 0) Heisenberg fixed 
point u* H = 6 e/(n + 8), we arrive at 



£ -2 9r 



(3.123) 



which demonstrates that the regimes as function of n 
for the existence of a nontrivial fixed point g* 2 become 
inverted as compared to the equilibrium case: g^f = 
2 (n — 4) e/n (n + 8) > only for n > 4, but is clearly 
unstable. The stable RG fixed point is thus character- 
ized by vanishing coupling <7q 2 = to the conserved field, 
which again implies decoupled model A dynamic critical 
behavior, with n = 0, v~ l = 2 — (n + 2) e/(n + 8), and 
zs = 2 to one- loop order, and purely diffusive z p = 2. 
For n < 4, on the other hand, gQ 2 — becomes unstable. 

Next, for = oo, the effective dynamic coupling in 
Eqs. 1|3.111[) - H3.113I) becomes 



9r = ®R9r ■ 



Thus Yd = and 



(3 W = 



w R9r 
(1 + wrY 



(3.124) 



(3.125) 



whence we see that w* D = oo is stable, which immediately 
implies that 7J = 7- = = 7^ = 7-, and 7 * = 7- = = 

7~ as well. Since 7* in Eq. 13.110(1 too reduces to the 
standard static expression, see Eq. I(3.75|l . this fixed point 
describes mere model A critical scaling, independent of 
the values of v* and g* 2 . In fact, 



Ps = 9r 



ur 



(3.126) 
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in addition to Eq. (|3.121[l , which only allows for the stan- 
dard decoupled model A Heisenberg fixed point u* H = 
6 e/(n + 8), g^ 2 = if n > 4. As for Q% = 0, there exists 
no finite nonequilibrium fixed point for n < 4. 



At the unstable fixed point 6* 



oo, Wq = 0, 



7s = 


0, 7x = 7~ 




7p = 


* * 


-- lb 




-2 


n + 2 


u* - 


+ 6 


h = 
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-e+ — 


+ 8 
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n 

-e+ — 


+ 2 
3 



u R -2g% 



(3.127) 
(3.128) 

(3.129) 
(3.130) 

(3.131) 



Again, the equilibrium relations l|3.72|) are satisfied. As 
to be expected, the decoupled model A fixed point u* H = 
6e/(n + 8), g^ 2 = is stable for n > 4 in the (u R ,g R ) 
subset of parameter space, whereas for n < 4 a novel 
fixed point u* x = 12e/(5n + 4), g*J = (4-n) e/2 (5n + 4) 
emerges, with unusual scaling exponents 77 = 0, 2/ = 
2-3 (n+4) e/2(5n + 4),z s = 2+ (4-n) e/2 (5n+4), and 
z p = 2 (all to one- loop order). But recall that this fixed 
point is unstable both in the wr and the Qr directions. 

This leaves us with the case Wq = 0, which according 
to Eq. I|3.114l) allows any value for the nonequilibrium 
parameter Qr. Yet since 



d(3 w _ 2 (n 
d^R =9 *\2- e « 



(3.132) 



at Wq — 0, stability requires Or < n/2. Quite generally, 
as in equilibrium, the model A fixed point (|3.94(l with 
decoupled diffusive relaxation for the conserved scalar 
density is stable for n > 4, with arbitrary 87? < n/2, 
but unstable for n < 4. The corresponding anomalous 
dimensions become 

7s = 0, ^ =7 i = - 7l = g**Q R , (3.133) 



7~=7 D =7 5 = 2 



r*2 



(3.134) 
&r ) ■ (3.135) 



7p 
T* = -2 



Remarkably, the equilibrium relations (|3.72() hold once 
again at any such fixed point. Inserting Wq — into 
Eqs. (|3.112|) and (|3.111|l and searching for nontrivial ze- 
ros leads to a quadratic equation, which is solved by 



9n 



e 



-3n(l -2Q R ) 



and ur 
with A : 



± ^9 n 2 (1 - 2 6t0 2 + 4 (71 - 4) 4 

"i + ^^ + e 



3e 
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n 2 (77 + 
16~ 



e \4 
+ (577 + 4) 9^(1-9^) 



(3.136) 
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FIG. 1: The nontrivial fixed points g 2 = gjj/e and u = un/e 
[Eq. (I3.136H as functions of Gh < n/2 for n = 2 (0.94 < 
9 fl < 1) and 7i = 3 (0.84 < Qr < 1.5). 



With appropriate sign choices, this reduces to the special 
cases with 9* = 0, 1, and 00 already explored above; e.g., 
for 9* q = 1, one finds g*J — 2 (4 — 71) e/n (71 + 8) and u* = 
24e/7i(ri + 8), i.e., u* = u* - 3 = u * H = 6e/(n + 8). 
Yet, as depicted in Fig. ^ Eq. (|3.13t)|l permits an en- 
tire interval of nontrivial fixed point solutions, namely 
for 0.94 < Qr < 1 for n = 2, and 0.84 < Qr < 1.5 
for 71 = 3. Apparently therefore, there exists a line of 
fixed points that describes slight perturbations from equi- 
librium for two- and three-component order parameters. 
Whereas the ensuing anomalous dimensions satisfy the 
equilibrium constraints (|3.72(l . the critical exponents in- 
dividually differ from their equilibrium model C values, 
and vary continuously as functions of the nonequilibrium 
parameter 9^, as shown in Fig. We cannot exclude, 
however, that this unusual feature might merely consti- 
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FIG. 2: Critical exponents u, zs, and z p for the isotropic 
nonequilibrium model C with n — 2 (0.94 < On < 1) and 
71 = 3 (0.84 < 0_r < 1.5) in the weak dynamic scaling regime 
with (wq — 0) as functions of Qr < n/2. 
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tute an artifact of the one-loop approximation. 

In summary, for the case of isotropic detailed bal- 
ance violation in model C with a scalar order parameter 
(ri = 1), no stable genuine nonequilibrium fixed points 
are found. The RG flow then takes the system to the equi- 
librium model C fixed point with strong dynamic scaling 
zs = z p = 2+a/v {vJq — 1), and the standard scaling ex- 
ponents as given in Sec. lIII CI However, for model C with 
two- or three-component order parameter, in equilibrium 



governed by weak dynamic scaling (zg < z p , 



0), at 



least to one-loop order lines of nonequilibrium model C 
fixed points are found that include the equilibrium fixed 
points, yet allow for continuously varying static and dy- 
namic critical exponents. For n > 4, the conserved scalar 
density effectively decouples from the order parameter, 
which then follows model A dynamic critical behavior. 
The effective noise temperature ratio naturally plays 
no role in this decoupled scenario. 



IV. ANISOTROPIC VIOLATION OF DETAILED 
BALANCE IN MODELS C AND D 

A spatially anisotropic nonequilibrium perturbation is 
applied to models C and D by dividing our ci-dimensional 
space into two sectors of dimensionality d\\ and d± (with 
dii + d± = d), and assigning to them different noise 

strenghts Dm and D±, respectively, for the conserved en- 
ergy density: D V 2 -> D\\ Vjj +D± Vj_, whence Eq. (|2~%|) 
is replaced with 

(7?(x,t)r ? (x',t'))= (4.1) 

-2 (d,, Vjj + Dj_ Vi) 5(x - x') S(t - t') . 

The conserved field noise in the two sectors can thus be 
thought of as being coupled to thermal reservoirs with 
different effective temperatures Tji and Tj_, whence we 
obtain two distinct nonequilibrium parameters 0ii/_l = 

Dun A/DA. Correspondingly, a new degree of freedom 
enters the problem in the form of 

o- = Q\\/e ± , (4.2) 

the ratio of the temperatures of the heat baths cou- 
pled to the conserved density in the two spatial sectors, 
each measured with respect to the order parameter (crit- 
ical) temperature. An anisotropic perturbation clearly 
requires that a ^ 1. We may choose the label assign- 
ments such that Oil < 0j_, i.e., < a < 1. In general, 
we must allow for the anisotropic noise in Eq. JUT} to 
induce further splittings in the renormalized parameters 
Dfl, and Ai? as well (see Ref. 01). For model D, 
we may in addition impose anisotropic strengths in the 
conserved order parameter noise, 

<r(x,i)C /3 (x',t'))= (4.3) 
-2 (A,, Vjj + Aj_ Vi) <S(x - x') 8{t - t') 5 af3 , 

whereupon 0||/j_ = ^ ||/_l A/D Aj|/j_. 



A. The anisotropic nonequilibrium model C 

1. Renormalization to one-loop order 

For model C (a = 0), we merely need to replace 
Dq 2 — > D||q 2 + D^q^. Therefore, the only modifica- 
tions to the previous perturbation expansions occur in 
those diagrams which contain an internal conserved field 
propagator GS. (q,w). Eqs. flOSt and (|^1T|) still hold, 

which implies 



z p = zz 1 = z~ 



(4.4) 



to all orders in perturbation theory. Removing the loga- 



rithmic divergence in T~^(q, 0) gives to first order 
Z D = l--g 



(4.5) 



as in Eq. i|3.40[) ■ The results of renormalizing those quan- 
tities that do not involve taking derivatives with respect 
to the external momenta are quite similar to the pre- 
vious isotropic results. The effects of the anisotropy in 
these cases is a mere replacement of the factors 1 — 
with 1 — (<i|| 0|| + d± Q±)/d in the expressions for the 
renormalization constants. For example, the fluctuation- 
induced T c shift becomes 



1 



(4.6) 
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r c /(l +w) + p 2 



After rewriting in terms of the true distance from the 
critical point r = r — r Cl subsequent multiplicative renor- 
malization of r~ s (0,0) and 8 uj T~ s (0,uj) \ u =o leads to 



{Z~Z S ) 1/2 Z X Z T = \- 



6 
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{Z-Zs) 1 ' 2 = 1 
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(4.7) 



(4.8) 



which are just the straightforward generalizations of 
Eqs. (|3~53jl and (233- 

Yet the resulting expressions become more complicated 
when the quantities to be renormalized involve deriva- 
tives with respect to the external momentum. Consider 
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[r/(l + w)+p 2 ] 3 



evaluating the integrals at the normalization point r = 
ji 2 in dimensional regularization then yields 



(Z~Z S ) 1/2 Z. 



An /j 



i__t e|| _^e ± V^~ 



(1+w) 3 d ~" d 

dx/u 3 2 (©|| - ©-l) A d fi 



3d (1 + w) 3 
Next, from (0,0), we get 



Z~ Z~ = 1 + — , 
s a ^ 1 + w V d 



"a ) e 



(4.10) 



(4.11) 



and likewise the renormalization of the three- and four- 
point functions results in expressions that can simply be 
obtained from Eqs. i|3.47[l - H3.46[l through the substitu- 
tion 6 -> (dp 6|| + d_L &±)/d. 



Upon identifying Z~ 



Z\,,, ± and factoring from 



these products of Z factors, we obtain at last 
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(4.17) 
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n + 8 A^~ £ 
z «n/i — 1 6 " e 
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1 + w 
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1 + w 
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, ... \ dx/n 9j - 9± 

in addition to Eqs. I|4.4() and Q4.5J1 . 



(4.18) 



AdH 



2. Anisotropic RG fixed points 

We first focus on the nonequilibrium anisotropy pa- 
rameter a. Yet because the anisotropic contributions 
to the renormalization constants Z~ and Zj, , (even 

A||/_L 11/^ 

when these are not chosen identical) and Z~ = Zr>„ , . , 

respectively, are equal at least to one-loop order, its RG 
/3 function reads 

Pc = °r(id, -713,1 + 7A„ -7X }] 

-75 x +7r>x -7Ax +7Xj) =0 ' ( 4 ' 19 ) 

leaving the fixed point a* undetermined. Indeed, consid- 
ering the heat bath ratios in the two sectors separately, 
we find (omitting the || / + subscripts on the renormal- 
ized couplings Wr and g^) 



'6,,/, - 



2w R g R Q\\/± R 



lw± (l + w R ) 2 

(4.20) 

in almost obvious generalization of Eq. (|3.114|l . and the 
corresponding stability matrix becomes 



<9/fo„/0||fi 3/3en /@±r 
dpeJQ\\ R d/3 e± /e^R 

d- 2d|| Q\\r - d±_ <d ±R 
-d|| 0±r 



x (4.21) 



d(l+w R ) 2 

-d± <d\\ R 
d|| 0||_r - 2d ± Q±r 



We restrict our investigation to the case WRg R > 
here, since for weak dynamic scaling with Wq = 0, we 
already found genuine nonequilibrium behavior even for 
isotropic detailed balance violation in Sec. lIIIEl We then 
find that the fixed points ©ij/xo = ® an< ^ ©|/Xoo = 00 
are unstable, whereas there is a stable line of fixed points 
given by 



dii o; + d ± e* ± = d , 



(4.22) 



which incorporates the equilibrium fixed point 9J„„ = 

I C H 

1 = 0^ cq . Its stability matrix eigenvalues are and 1. 
The marginal flow direction is clearly along the fixed line 
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(arbitrary value a*). In fact, one readily computes for 
the anomalous dimensions 
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Consequently, the relations (|3.72() that follow from the 
fluctuation-dissipation theorems for the order parame- 
ter and conserved density, respectively, are fulfilled even 
here. In this sense, the entire fixed line (|4.22|) again rep- 
resents a system mimicking thermal equilibrium, albeit 
with potentially anomalous scaling exponents. 
Next we consider 
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Thus, the weak scaling fixed point Wq = is stable for 
n > 2 ± (2d^/||/3d) (6|| — 6j_), whereas w* D = oo is 
unstable in the w direction. In addition, there appears 
a strong scaling nontrivial fixed point w* given by the 
solution of the nonlinear equation 



_(l + tW .)--(l + ^ = ± - 3d 

Since at this fixed point 
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it is stable for n > 4/3 (1 + w*). The associated anoma- 
lous dimensions become 
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FIG. 3: The nontrivial fixed points g 2 = gjj/e and u — un/e 
[Eq. H4.40I 1 as functions of c g , 3n/4 < c g < l + 3n/4 for n = l, 
2, and 3. 



We now need to find the RG fixed points from the coupled 
(3 functions 
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Note that 3n/4 < c g < 1 + 3n/4. As usual, the model 
A fixed point u* H = 6e/(n + 8), #q 2 — is unstable for 
n < 4, but becomes stable for n > 4. For n < 4, the 
stable RG fixed point acquires nonzero values for both 
couplings ur and 5^, to be found as the solution of the 
quadratic equation following from Eqs. (|4.36|) and H4.37|) : 
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e 

Jc 

3e 



B ± ^B 2 - 8 (4 - n) C 
ZcgffR 
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(4.40) 



9 • 



n + 2 \ e 
with B = (n + 2) (8 - n) - 4 (4 - n) c 
and C = 2 (n + 2) 2 - (n + 2) (8 - n) c g + 2 (4 - n) c 2 . 

This incorporates the equilibrium fixed point, since with 
w* c — (2/n) — l one obtains B — 6n and C — n 2 (n+8)/8, 
whence the two solutions in Eq. H4.40J1 reduce to the two 
nontrivial model C fixed points discussed in Sec. IIII CI 
gf = 2e/n, u* = and g*J = 2 (4 - n)e/n{n + 8), 
u* H = 6 e/ (n+8) for g 2 ^ and ur — ur—3 g R . The solutions 
(|4.40|) for n = 1,2,3 are shown in Fig.|31as functions of c g . 
Upon inserting into the anomalous dimensions (|4.33|) - 
(|4.35|l , this once again yields continuously varying static 
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FIG. 4: Critical exponent n for the anisotropic nonequilib- 
rium model C with n = 1, 2, 3 as functions of the parameter 
c g (3n/4 < c 9 < l + 3n/4). 



B. The anisotropic non-equilibrium model D 

We next consider the critical behavior of our nonequi- 
librium version of model D with spatially anisotropic 
conserved noise. The anisotropy may now be imposed 
through Aq 2 — > Ay q| + Aj_ q^ in addition to Dq 2 — > 

D|| q| + Dj^ qj_, see Eqs. <|4.3[) and (|4.1|) . We first com- 
pute the fluctuation-induced T c shift from the criticality 
condition 9 g 2 T~ s (q, 0) | 9 =o= 0. Yet for this purely re- 
laxational dynamics, at least to one-loop order the or- 
der parameter noise strengths Aim do not enter, and we 
again arrive at Eq. (|4.6|) as for model C with noncon- 
served order parameter. However, since asymptotically 
w — > oo here, the nonequilibrium parameters 6y/j_ dis- 
appear, simply leaving the static one-loop T c shift 



1 



(4.41) 



and dynamic critical exponents 77 = —75, v~ x = —7*, 
and zs = 2 + 7* = z p as depicted in Figs. 0] and \5\ as 
functions of the parameter c g . Via Eq. (|4.38() , the param- 
eter here is the time scale ratio w* , or equivalently, the 
effective temperature difference Oy — 0j_ between the lon- 
gitudinal and transverse sectors, see Eq. 14.31|l . In con- 
clusion, the one-loop RG flow equations for the nonequi- 
librium model C with spatially anisotropic noise allow for 
novel strong dynamic scaling regimes with z$ = z p with 
continuously varying critical exponents even for a scalar 
order parameter that encompass the equilibrium model 
C fixed point. 



In the same manner, all other fluctuation contributions 
reduce to the equilibrium expressions. As demonstrated 
explicitly in Sec. IjHI All by integrating out the conserved 
scalar density, the terms violating detailed balance be- 
come obsolete at the model D fixed point w* D = 00. We 
remark that generalizations of dynamical models with 
conserved order parameter that contain reversible mode 
couplings to other conserved quantities behave markedly 
different when subject to spatially anisotropic noise cor- 
relations: In models H and J with 'dynamical' noise, the 
nonlinear mode couplings induce anisotropic shifts of the 
critical temperature already to one-loop order, thus ren- 
dering the fluctuations soft only in one subsector of mo- 
mentum space. For the ensuing two-temperature models 
H and J, no stable RG fixed points could be identified, 
perhaps indicating that no simple nonequilibrium steady 
state is approached in the long-time limit 0, [l8L fli^ . 
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FIG. 5: Critical exponents v and zs 
nonequilibrium model C with n 

(3n/4 < c g < 1 + 3n/4). 
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for the anisotropic 
1,2,3 as functions of c„ 



C. The two-temperature model D 

1. Derivation of the effective theory 

The anisotropic non-equilibrium model D discussed 
in the previous subsection does not actually represent 
the most general spatially anisotropic extension of re- 
laxational dynamics with a conserved order parameter 
coupled to a conserved scalar density. Rather, one can 
generalize Eqs. I|2.5|l and (|2.6|l with a = 2 to 



dS a 



(4.42) 



A,, V 2 



11 L 

A,V 2 : 



r. - VI 



6 



E^ 2 



g\\ p 



s a 



s a + c , 
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dt 11 11 



P- 



(4.43) 



with noise correlators 

<r(x,t)c /J (x',o) = 



-2 (Ay Vjj + A ± Vi) S(x - x') <5(< - t') S ap 

(r,(x,<)77(x',<')> = 

-2 (5|| Vfj + 5_l Vi) S(x - x') <5(t - t' 



(4.44) 

3/3 

) 

(4.45) 

We choose the labels such that r± < rii, which is to 
be interpreted as a lower order parameter temperature 
j"-L < yll j n the transverse spatial sector. Thus, at the 
critical point, the longitudinal fluctuations remain un- 
critical ('stiff'), similar to equilibrium anisotropic elas- 
tic phase transitions [5l|, or the behavior at Lifshitz 
points [23. Nonlinearities and higher-order gradient 
terms should then only be relevant in the 'soft' trans- 
verse directions. In analogy with the two-temperature 
nonequilibrium model B (or randomly driven lattice 
gases) 0, [2 El ! it is possible to construct an effective 
field theory which reduces our most general anisotropic 
model D to an equivalent equilibrium system, albeit with 
spatially long-range correlations. We first construct this 
effective field theory and then perform the p erturba- 
tional renormalization of the model to one-loop order, 
discussing finally the ensuing RG flow equations. 

Since T11 = r|| — r c \\ > in the noncritical momen- 
tum space sector, whereas t± = rj_ — r c ± — ► at the 
phase transition, we expect the terms oc q|, q| to 

be irrelevant as compared to q^. In fact, in the Gaus- 
sian theory at criticality A|| Ty q 2 ~ q^. Hence we apply 
anisotropic scaling with [qj_] = p, [qy] = [q±] 2 = p? , 
[ui] = [q_i_] 4 = /i 4 , which yields the following scaling di- 
mensions: [Aj_] = [AjJ = p°, [A||] = [A||] = p~ 2 , A = p~ 4 , 

[ry±] = p\ [Dx] = [Dx] = M 2 , = [£>„] - pP, 
and with [S a ] = f r 1+d w+ d ^/ 2 , [p] = [/\\+ d ±/ 2 at last 
[uyx] — [flf/J = p 4 ~ d ~ dn - Consequently, the longitudi- 
nal parameters all become irrelevant under scale trans- 
formations, except the marginal product [Ay tii] = p . 
Therefore in the vicinity of the critical point, all nonlin- 
earities in the longitudinal sector and fluctuations ~ q|, 

q| qj_ can be safely omitted. From the naive scaling di- 
mensions of u\\/x an d 9\\/x we infer the upper critical 
dimension 

d c =A-d\\. (4.46) 

It is reduced as compared to the isotropic case be- 
cause the critical fluctuations are confined to the d±- 
dimcnsional subsector here. 

To proceed further, we rescale the fields according to 



S a -» (X ± /X ± )^ 2 S a , p -» (A^/A^) 1 / 2 ^ and define 



A 



A 



A 



(4.47) 



The effective Langevin equations of motion for the order 
parameter fields S a and the conserved density p near the 
phase transition at last read 



dS a 



= A; 



cV 2 +Vi (rx-Vi) 



,S' a 



-A, V 2 , 



ux 
6 



% - ^ 



gxp 



(4.48) 

s a + c , 

(4.49) 



with the corresponding noise correlations 

(C(x,i)C /3 (x',<')> = 

-2 A ± Vi (5(x - x') 5{t - t') 5 af} 

(7 7 (x,t)7 ? (x , ) t')) = 

-2 D ± 6x Vi 5(x - x') S(t - t') 



(4.50) 



(4.51) 



where Q± again denotes the heat bath temperature ratio, 
9-L = ^ ^ • (4-52) 



Dx Aj 
Dx A 1 



The preceding Eqs. H4.48fl - H4.51fl define the two-tempera- 
ture nonequilibrium model D. Our analysis that led to this 
effective critical theory for the most general nonequilib- 
rium model D with dynamical anisotrop y cl osely parallels 
that of the two-temperature model B |n|, 0, 131. No- 
tice that after the field rescaling, only the noise strength 
in Eq. (|4.51|) violates the Einstein relation with the cor- 
responding relaxation constant Dx in the critical trans- 
verse sector, if Ox 7^ 1. 

Yet we can certainly write Eqs. I|4.48[) and (|4.49|) in the 
form of purely relaxational Langevin dynamics 



8S a (x,t) 
dt 
<9p(x, t) 



x y2 5n eS [S,p} 



\x vj_ 
DxVl 



6S a (x, t) 
2 5H c fi[S,p] 



dt ^ ^ 6p(x,t) 

with an effective Hamiltonian 



+ C(x,i) , (4.53) 
+ *7(x,t), (4.54) 



H e ff[S,p] 

d d q 



(27T) 

d d x 



d 



qi(n 



ux 
4! 

91 



2ql 

£S Q (x) 2 ^(x) 

a, f3 



S a (q)S a (-q) 



p(x)^^(x) 2 + -p(x) 



(4.55) 
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As is obvious from its harmonic part, this Hamiltonian 
contains long-range interactions generated by the dy- 
namical anisotro py, akin to those found in driven dif- 
fusive systems El llll. but also at equilibrium ferroelastic 
phase transitions [21j and at Lifshitz points |22|. How- 
ever, our earlier investigations of model D subject to var- 
ious nonequilibrium perturbations showed that the heat 
bath temperature ratio H4.52|) disappeared entirely in the 
asymptotic limit; since we are here concerned only with 
the transverse sector, this is reached as wj_ = D^/X^ — > 
oo. Indeed, integrating out the conserved scalar den- 
sity from the dynamical action proceeds precisely as in 
Sec. IIII Al since the dynamically generated long-range in- 
teractions only appear in the order parameter propaga- 
tor. As a result, the remaining detailed balance violation 
plays no role at all for the fixed point properties, and the 
two-temperature model D in effect becomes an equilib- 
rium system. Thus, we expect it to relax towards a sta- 
tionary state that is characterized by the Gibbsian prob- 
ability distribution "P C q[S, p] cx exp(— H c s [S, p\), with the 
effective long-range anisotropic Hamiltonian l|4.55|l . 



2. Renormalization and critical exponents 

We introduce the renormalized fields and parameters 
as in Sec. IIII Bl supplemented with 

c a = Z c c , Q± R = Z e± 6_l • (4.56) 

But the deviation from the critical dimension now reads 



= d r - d = 4 - d - d\\ = 4 - 2 d\\ - di 



In summary, we obtain the exact relations 



Zs = Z- 1 = Z\± — Z 1 



Z n = Z. 



s 
-l 



Z 



(4.60) 
(4.61) 



The perturbation expansion naturally acquires the 
same structure as for the equilibrium model D (or model 
C, with w — > oo). To one-loop order, which is determined 
entirely by simple combinatorics, we can in fact immedi- 
ately take over the equilibrium renormalization constants 
(|3.62l) - (|3.67[) with shifted critical dimension, the replace- 
ments u — > u±, g 2 — > (jrj_/c d ii/ 2 , an d modified geometry 
factor Ad — > A(d\\,d±) as given in Eq. I|4.58|l . This is 
confirmed explicitly by renormalizing 



r~ s ( q ,o) = A 
x(u x - 3gi) 



cqf +ql 



/ , n + 2 

q± r± i 



6 



Pi 



(4.62) 



/plcpf + Pi^ + p^) 2 ] 2 

at the normalization point tj_ = fi 2 , which leads to 

Z\ ± = 1 j (4.63) 
n + 2 u^-3g\ A(d\\,d±)n~ 



Zt 1 6 C <V 2 e 
Similarly, the logarithmic singularity in 



(4.64) 



r~(q,0) = D ± qj 



pi 



P [cpf + Pl(r±+p±) 2 ] 2 



(4.57) 

is absorbed into 



(4.65) 



and we define the anisotropic geometric factor as 



A{d\\,d_[ 



r(3-d/2-d||/2)r(d/2) 

2 rf - 1 7T d / 2 r(d J ./2) 



(4.58) 



with A(Q,d) = Ad- As a consequence of the conserva- 
tion laws, and the ensuing momentum dependence of 
the vertices, in analogy with the isotropic situation (see 
Sec. IIII B|) the following relations hold to all orders in per- 



turbation theory: 
-2D ± 9j , r~ 



and d n 2 r~~(qj_,0) 
1, Z 2 



q±=a= 

l<?J_=0 = 

Zs=h 



r~ (o,w) 
pp v ' 1 

ss (0,u) = iw, — - r± - ss 
-2 A ± , whence Z~Z p = 1, Z~Z D± Z e± = ^ ^ s 
and Z~Z\ ± = 1. Notice that since the order parameter 
Langevin equation fulfils the Einstein relation, this sat- 
isfies the identity Z\ ± = {Z$ / Z-) 1 / 2 following from the 
fluctuation-dissipation theorem 



(q, uj) 



2A 



Imr~ g (q,w) 



(4.59) 



Moreover, none of the nonlinear vertices carries trans- 
verse momentum, which leaves the c q| term in the prop- 
agator unrenormalized to all orders in perturbation the- 
ory as well, d q 2 r~ s (q||,0) | g|l =o= A_l c, i.e., Z x± Z c = l. 



n g\ A{d h d 1 _)ix' 
Zd ± = i — — 



(4.66) 



2 C <V 2 e 
Finally, the three- and four-point vertex functions yield 
n + 2 u± A(d\\ , d±) 



Z g = l- 



Z u = l- 



3 c d \\l 2 e 
n + 2 g\ A(d\\,d±)fi- 



2 c d ii/ 2 e 
n + 8 u± A(d\\,d±) fi~ 



6 



od||/2 



6 1- 



g±\ 9± Md\\,d±)(jr 



(4.67) 



(4.68) 



as well as Z p = Zp, whence Zq ± = 1 as expected: Since 
the nonequilibrium parameter Q± disappears from the 
asymptotic theory entirely, its fixed point remains unde- 
termined, /3q ± = 0. We also remark that the T c shift 
obtained from the criticality condition is 



Foe 



(n + 2) (uj. -3gj) A(d||,d^. 
3c d ii /2 (d + d\\ -2)(4-d-d||) 



1 2/(4— d-dn) 



(4.69) 
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The divergence of the denominator here indicates that in 
addition to the reduction of the upper critical dimension, 
the lower critical dimension is lowered as well to d\ r = 
2— du , just as in the two-temperature model B [lttliull2| . 

As in Sec. IIII B 31 we can now define flow functions via 
logarithmic derivatives of the Z factors with respect to 
the renormalization scale /i, see Eas. (|3.57(l - <|3.59[1 . with 
{a} = u±, g±, Aj_, c, and t± here. The solutions to 
the RG equations for the vertex functions are given by 
Eq. I|3.61[) . with running couplings and parameters deter- 
mined by the flow equations Eq. (|3.60|) . The general scal- 
ing form for the renormalized order parameter response 
and correlation function thus obtained at an IR-stable 
fixed point becomes 



x(Tx,qn,q_L,w) = q ± 2+n x 



C(t±, q||,q_L,w) = q ± 



q± 'l. 

_ -2—z+T) At r J- 



1± — ) ( 4 70) 

1/!/' „1+A' z I' \' i - nJ ) 



'C 



l/u ' 1+A ' „z 



(4.71) 



where, in addition to the usual static exponents 77, v, and 
the dynamic exponent z, the anisotropy scaling exponent 
A has been introduced. 

Since the two-temperature model D is effectively in 
equilibrium, we may insert the Heisenberg fixed point 
u* H = 6 e/ [n + 8) for u± = u± — ig\ to obtain the static 
critical exponents, which thus assume the usual one-loop 
form 



V = -7s = 



= -7;, = 2 - 
2 - dv 



n 



" n + 8 

4 - n 



2(n + 8) 



(4.72) 
(4.73) 

(4.74) 



but with e = 4 — d — d» . To two-loop or der, the 
static critical exponents were evaluated in Ref. [13j. In 
addition, upon invoking the exact relation l|4.60|l . i.e., 

7s = J\ ± = —7c w e arrive at 



z p = 2- 
zs = 4 



= 4-77 



A = l--i£- = l-^ 
2 2 



(4.75) 
(4.76) 

(4.77) 



All order parameter scaling exponents are thus given 
by the static critical exponents, precisely as in the two- 
temperature model B [ill H2I Il3l | . 

As in equilibrium, the same is true for the dnyamic 
critical exponent governing the conserved energy density 
p. Since £ c = in the one-loop approximation, the en- 
suing RG (3 functions for u± and are just Eqs. (|3.78|) 
and H3.79H of the equilibrium model C/D. Consequently 
for n < 4, g\ -> g* c 2 = 2 (4 - n) e/n(n + 8) = 2a/nv 
and 



Q 
V 



(4.78) 



whereas for n > 4, a < and lj\ — > 0. Therefore the 
coupling between the order parameter and the conserved 
density becomes irrelevant, resulting in a purely diffusive 

z p ee 2 . (4.79) 

Therefore, the independent static and dynamic critical 
exponents to one-loop order look identical with those of 
the equilibrium model D, albeit with shifted e = 4 — d — 
d\\ . The order parameter scaling exponents, including the 
additional anisotropy exponent, moreover are precisely 
those of the two-temperature model B. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have studied the critical behavior 
of the relaxational models C and D with nonconserved 
and conserved order parameter, respectively, coupled to 
a conserved scalar density, and subject to both isotropic 
and anisotropic nonequilibrium perturbations. This sup- 
plements previous work on the identification of genuine 
nonequilibrium critical behavior in the form of modified 
dynamic universality classes in 0(n)-symmetric models 
|19i | . These investigations have demonstrated the general 
robustness of the equilibrium critical behavior in mod- 
els with nonconserved order parameter with respect to 
the violation of detailed balance, both isotropically and 
anisotropically. This remarkable stability has been es- 
tablished particularly for model A which represents the 
simplest critical dynamics with a nonconserved order pa- 
rameter 0, H, 13 But even in more complicated sit- 
uations involving reversible mode couplings between a 
nonconserved order parameter and additional conserved 
quantities, viz. models E, G or their n-component gen- 
eralization, the SSS model |2^|, the equilibrium RG fixed 
point turned out to be stable, and thus describes the 
asymptotic critical power laws, despite the existence of 
additional genuine nonequilibrium fixed points [lflj |. 

Our results here for model C with scalar order pa- 
rameter (n = 1), which extends model A to include a 
nonlinear coupling to a conserved scalar density, are in 
accord with this general observation. Specifically in the 
case of isotropic detailed balance violation, the coupling 
of the order parameter and the conserved field to differ- 
ent heat baths gives rise to the nonequilibrium param- 
eter which represents the temperature ratio of these 
heat baths. This variable induces different renormaliza- 
tions for the noise strenghts, with the possibility for gen- 
uinely new dynamic as well as static critical behavior. 
Even when unstable, such nonequilibrium fixed points 
would affect crossover features and corrections to scaling 
in the critical regime. However, a stability analysis yields 
that the equilibrium fixed point that describes strong dy- 
namic scaling (wq = 1) with = 1 remains stable. At 
least to one-loop order, we could not identify any gen- 
uine nonequilibrium model C scaling regime for the case 
of a scalar (Ising) order parameter, even for the extreme 
situations with either = or = 00. For n = 1, 
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the asymptotic critical behavior is thus definitely gov- 
erned by the equilibrium model C fixed point, with the 
static critical exponents of the 0(n)-symmetric </> 4 model, 
and with equal dynamic exponents z$ — z p = 2 + ot/v 
[TH Il6j . The critical behavior again reduces to that of 
the isotropic case as described above. Therefore, we ob- 
tain the remarkable result that a quadratic coupling of a 
scalar order parameter to a conserved density, which pre- 
serves the internal symmetry of the corresponding equi- 
librium system, does not produce any novel universality 
classes for models with a nonconserved order parameter, 
subject to detailed balance violations. This result is to 
be seen in contrast with the system studied in Ref. [24| 
which incorpates a linear coupling of a conserved field to 
a nonconserved order parameter; in that case, effective 
long-range interactions are generated, which yield novel 
nonequilibrium scaling features. 

For model C with n-component order parameter, RG 
fixed points with 0* ^ 1 do appear for n > 4. Yet in 
this situation, the order parameter effectively decouples 
from the conserved density, resulting in model A criti- 
cal behavior. However, our RG analysis yields more in- 
teresting critical scaling for the nonequilbrium model C 
with two or three order parameter components. In equi- 
librium, one encounters weak dynamic scaling in these 
cases, with z$ < z p (wq = 0). We find that the one- 
loop flow equations allow for an entire line of fixed points 
encompassing the equilibrium case. Consequently, an 
interval of fixed point values emerges for the nonequi- 
librium parameter 0, leading to continuously varying 
static as well as dynamic critical exponents (as shown in 
Fig. [2J • Curiously, the general scaling relations imposed 
by the fluctuation-dissipation theorem remain satisfied 
along this entire fixed line. In a similar manner, for the 
nonequilibrium model C with spatially anisotropic, dy- 
namical noise, we obtain a line of nonequilibrium strong 
scaling fixed points for n < 4, i.e., even for a scalar order 
parameter, with an allowed interval of fixed point values 
w* , characterized again by continuously varying scaling 
exponents (Figs. 0] and [SJ. 

For the nonequilibrium model D with isotropic detailed 
balance violation, the relaxation of the conserved non- 
critical density occurs inevitably much faster than that 
of the also conserved order parameter. Hence the con- 



served energy density is always able to keep up with the 
critical fluctuations and does not in turn influence the 
order parameter dynamics. This is clearly seen in the 
perturbation expansion upon taking the limit w — > oo 
for the diffusion rate ratio of the conserved scalar den- 
sity and order parameter, whereupon all terms involving 
the heat bath temperature ratio disappear. A further 
rescaling of the static nonlinearity u — > u and the cou- 
pling constant g 2 — > g 2 then reduces this nonequilibrium 
model D variant fully to its equilibrium counterpart. 

However, introducing dynamical anisotropy, i.e., dif- 
ferent effective noise and ordering temperatures in the 
longitudinal and transverse spatial directions in model D 
with conserved order parameter has a much more drastic 
effect, since now only the momentum space sector with 
weaker noise softens. As with the anisotropic nonequi- 
librium model B 0, 0, ^| , it is possible to recast the 
emerging two-temperature model D with its nonlinear 
coupling to a conserved density into an effectively equi- 
librium model, albeit with a Hamiltonian that already 
contains long-range correlations. The consequences are 
strongly anisotropic scaling, and a reduced upper critical 
dimension d c = i—dn . We finally remark that this feature 
of the two-temperature relaxational models B and D is at 
variance with other conserved order parameter systems 
that incorporate reversible mode couplings to additional 
slow variables. Upon introducing anisotropic dynamical 
noise into models J [13 and H the equilibr mm m- 
tegrability conditions become irretrievably violated; at 
least to one-loop order one cannot even find any stable 
RG fixed points, suggesting that simple nonequilibrium 
steady states may not be sustainable in those situations. 
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APPENDIX A: EXPLICIT ONE-LOOP RESULTS FOR THE VERTEX FUNCTIONS 



In this appendix, we list the explicit results to one-loop order in the perturbation expansion for the vertex functions 
required for the renormalization of the parameters and couplings of models C (a = 0) and D (a = 2). In the following 
expressions the momentum integrals are given in abbreviated notation, i.e., f . . . = (27r) _d Jd d p . . ., and the internal 
frequency integrals have already been performed (via the residue theorem). We do not provide the Feynman graphs 
themselves, but only note the number of the contributing one-loop diagrams for each vertex function. 

For r~ s (q, uj), there are three one-loop graphs that give 

r ?o(q^) = Aq' 



n + 2 ,~ „~o, f 1 2 
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Only one one-loop diagram contributes to each of the other three two-point functions. The resulting expressions read 
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There arc three one- loop diagrams that contribute to the three-point function r~ ss (q, u>). Here, q and w denote 
the wavevector and frequency of the outgoing external p leg. The vertex function is evaluated at symmetric incoming 
labels — q/2 and — ui/2 for the order parameter fields. Setting the external frequency cj to zero, we obtain 



r~ ss (q,0) = Dq 2 S 
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In the same notation, four one-loop graphs yield 
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Finally, we need the four-point vertex function r~ ggg (q, 0), for which there are ten one- loop Feynman diagrams. 
We merely record the final result for q — > 0; after a little algebra, one arrives at 
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+3g 2 <3 



p p a (r + p 2 ) + f p 2 \r + p 2 p a (r + p 2 ) + f p 2 



+3g 2 



A 
^2 



A i P (r + P 2 )[p"(r + P 2 ) + jP 2 ] \r + p 2 p" (r + p 2 ) 



2\ i D n2 



-sSe 



?p 2 



u Jp p a (r + p 2 ) + f p 2 V r + p 2 p a (r + p 2 ) + f p 2 (r + p 2 ) [p a (r + p 2 ) + f p 2 ] 



(A7) 



5 4 1? 



#P 2 



" A ip(f + p 2 )[p a (r + p 2 ) + f p 2 ] U+p 2 p a( r + p 2) + D p 2 ( r + p 2)[pa( r + p2) + C p2 ] 



[1] P.C. Hohenberg and B.I. Halperin, Rev. Mod. Phys. 49, 
435 (1977). 

[2] R. Graham, Springer Tracts in Modern Physics, Vol. 66 
(Springer- Verlag, Berlin, 1973); U. Deker and F. Haake, 
Phys. Rev. A 11, 2043 (1975). 

[3] Prominent counter-examples emerge in certain reaction- 
diffusion systems, see, e.g., U.C. Tauber, Acta Phys- 
ica Slovaca 52, 505 (2002); U.C. Tauber, e-print 
cond-mat/0304065 (2003), and references therein. 

[4] H.K. Janssen and B. Schmittmann, Z. Phys. B 64, 
503 (1986); K.-t. Leung and J.L. Cardy, J. Stat. Phys. 
44, 567 (1986); for a comprehensive overview, see: 

B. Schmittmann and R.K.P. Zia, in: Phase Transi- 
tions and Critical Phenomena, edited by C. Domb and 
J.L. Lebowitz, (Academic Press, London, 1995), Vol. 17. 

[5] M. Kardar, G. Parisi, and Y.-C. Zhang, Phys. Rev. Lett. 
56, 889 (1986); for recent reviews, see: J. Krug and 
H. Spohn, in: Solids far from Equilibrium: Growth, Mor- 
phology and Defects, ed. C. Godreche (Cambridge Uni- 
versity Press, Cambridge, 1990); T. Halpin-Healy and Y.- 

C. Zhang, Phys. Rep. 254, 215 (1995); A.L. Barabasi and 
H.E. Stanley, Fractal Concepts in Surface Growth (Cam- 
bridge University Press, Cambridge, 1995); J. Krug, Adv. 
Phys. 46, 139 (1997). 

[6] This is the case, e.g., for the anisotropic KPZ equation: 

U.C. Tauber and E. Frey, Europhys. Lett. 59, 655 (2002). 
[7] F. Haake, M. Lewenstein, and M. Wilkens, Z. Phys. B 

55, 211 (1984). 
[8] G. Grinstein, C. Jayaprakash, and Y. He, Phys. Rev. 

Lett. 55, 2527 (1985). 
[9] K.E. Bassler and B. Schmittmann, Phys. Rev. Lett. 73, 

3343 (1994). 

[10] U.C. Tauber and Z. Racz, Phys. Rev. E 55, 4120 (1997); 

U.C. Tauber, J.E. Santos, and Z. Racz, Eur. Phys. J. B 

7, 309 (1999); e: ibid. 9, 56 (1999). 
[11] B. Schmittmann and R.K.P. Zia, Phys. Rev. Lett. 66, 



357 (1991). 

[12] K.E. Bassler and Z. Racz, Phys. Rev. Lett. 73, 1320 

(1994); Phys. Rev. E 52, R9 (1995). 
[13] B. Schmittmann, Europhys. Lett. 24, 109 (1993). 
[14] H.K. Janssen, Z. Phys. B 23, 377 (1976); C. De Domini- 

cis, J. Phys. (Paris) Colloq. C 1, 247 (1976); R. Bausch, 

H.K. Janssen, and H. Wagner, Z. Phys. B 24, 113 (1976). 
[15] B.I. Halperin, P.C. Hohenberg, and S.-k. Ma, Phys. Rev. 

B 10, 139 (1974); ibid. 13, 4119 (1976). 
[16] E. Brezin and C. De Dominicis, Phys. Rev. B 12, 4954 

(1975). 

[17] R. Folk and G. Moser, Phys. Rev. Lett. 91, 030601 
(2003). 

[18] J.E. Santos and U.C. Tauber, Eur. Phys. J. B 28, 423 
(2002). 

[19] An overview is given in: U.C. Tauber, V.K. Akkineni, 
and J.E. Santos, Phys. Rev. Lett. 88, 045702 (2002). 

[20] P. Sen and S.M. Bhattacharjee, J. Phys. A 35, L141 
(2002) 

[21] R.A. Cowley, Phys. Rev. B 13, 4877 (1976); R. Folk, 
H. Iro, and F. Schwabl, Phys. Lett. 57 A, 112 (1976) 
112; Z. Phys. B 25, 69 (1976). 

[22] R.M. Hornreich, J. Magn. Magn. Mat. 15-18, 387 (1980); 
W. Selke, in: Phase Transitions and Critical Phenom- 
ena, edited by C. Domb and J.L. Lebowitz, (Academic 
Press, London, 1992), Vol. 15; for recent developments, 
see: H.W. Diehl and M. Shpot, Phys. Rev. B 62, 12 338 
(2000); H.W. Diehl, Nucl. Phys. B 612, 340 (2001). 

[23] L. Sasvari. F. Schwabl, and P. Szepfalusy, Physica 81 
A, 108 (1975); H.K. Janssen, Z. Phys. B 26, 187 (1977); 
C. De Dominicis and L. Peliti, Phys. Rev. B 18, 353 
(1978). 

[24] G. Grinstein, C. Jayaprakash, and J.E.S. Socolar, Phys. 
Rev. E 48, R643 (1993). 



